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Abstract 

We estimate possible systematic error that can arise in the simph- 
fied simulation of the Wr production and its subsequent decay into heavy 
neutrino at LHC (left- right symmetric model). The general- purpose event 
generators simulate this process as a Wr production followed by two de- 
cays, loosing full information about the polarization of the intermediate 
particles. We calculate the cross section using the full matrix elements 
with propagators, perform the simulation and compare the results with 
the simplified simulation. 
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1 Introduction 

It is known that in the Standard Model there are left-handed W-bosons 
that interact with left-handed fermions. The Left-right symmetric model 
(for more details see jT]) states that there are also right-handed W bosons 
interacting with right handed-fermions. In the minimal Left-right sym- 
metric model they have the same coupling constants to the fermions 
but are much heavyer, otherwise they would be already discovered. In 
other words, in the SM the terms responsible for the interaction between 
fermions and weak bosons are 



Here Q are quark doublets and i/) are lepton doublets. The Left-right 
symmetric model contains also the following terms: 



The lepton doublets here contain ordinary charged leptons and heavy 
right-handed neutrinos. We assume that they can be Dirac or Majorana 
neutrinos. 

The searches for these new particles are being performed now at LHC 
[2]. The purpose of this work is to estimate possible systematic errors 
arising from some approximations usually used in the simulation of the 
signal by general purpose event generators. 

We consider the following process, d-quark and u-antiquark interact 
with each other producing the -boson. It decays into the lepton LI 
and heavy antineutrino A*' that decays into another lepton L2, it-antiquark 
and d-quark via the virtual W^-boson. This process is represented in 
Figure [T] 

This process takes place regardless of the type of the neutrino: it 
can be either Dirac-like or Majorana-like. However if the neutrino is 
Majorana-like, another process is possible with the same initial particles: 
du L\L2ud (it is drawn in Figure [21). 

The cross-section of these reactions can be calculated explicitly. How- 
ever, the 4 - particle final state is too complicated for the general purpose 
event generators (for example pythia6 0) usually used in the HEP analy- 
ses. In such generators this process is considered in the following approx- 
imate way: quark-antiquark pair produces real W^~-boson, then it decays 
into real lepton L2 and antineutrino and the antineutrino decays into an- 
tilepton L2 and quark-antiquark pair (if the neutrino is Majorana-like it 
can decay either to lepton or to antilepton). In this case the information 
about the polarization of intermediate particles is lost. We consider both 
approaches and compare the results they give. Let us note that we neglect 
other possible reactions with the same initial and final state particles. 
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2 Cross section with the full matrix ele- 
ment 

Let the particles in our process have the following momenta: 

pi -momentum of the initial d-quark 

P2 -momentum of the initial it-antiquark 

ki -momentum of the lepton LI 

k2 -momentum of the lepton L2 

qi -momentum of the final d-quark 

52 -momentum of the final it-antiquark 

The amplitude of the reaction represented in Figure [T] is written in the 
following way (here we use the common notations that can be found for 
example in [3], [5]): 

ig ^(1 + 7^) / ^ ^* f (Pi +P2)m(pi +P2) 



^(,, N W -,^ (1 + 7') »(pi+P2-fci)+m, ig ^^ (1 + ^5) 

{qi + q2)\{qi + q2)p\ ig . p(l + 7^) , x /g\ 
ffV T72 -i=<iin 7, "(92) (3) 



(gi + q2? - V J V2 

Here we do not take into account any possible mixings among the particles. 
Further we also neglect the masses of all particles except VK~ boson and 
the neutrino. Using the standard rules (see [1]) one can show that 

|..|2 ^ 4/(fcip2)(g2fc2)((fcip2)(pigi) + (pi(fci -p2))(gi(fci -P2))) . . 

' ■ ' ((Pl +P2)2 - M2,)2((gi + g2)2 - Af2,,)((Pl +P2 " ^l)' -ml) ^ ' 

We used the mathematica package feyncalc to obtain this expression. 

Then this formula is substituted into the general expression for cross- 
section of massless particles [6] 

J 1 xWfi. , ; , , MA^|2 1 d^ki<fk2d^qi<fq2 

da^j^S^\k, + k2+q^ + q2~Pi^P2)\M\ 4(p^p^) i6fcOfcO,o,o (5) 

The delta- function of spatial momenta is taken off by integration over q2. 
To take off the delta-function over energy we rewrite 

d'^qi = [qii^ d\qi\ sin OdOdif 

where 9 and (p is polar and axial angle relative to the vector A:i+fc2— pi— P2. 
We also rewrite the argument of the delta-function as 

\qi \ + \J |(?ip + |fci + — Pl — P2p + 2|(ji||fci + fc2 — Pl — P2I COS& + 

|fci| + Ifel — |pi| — |p2 
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After integration over \qi\ we obtain 

, ^ 4ff^(fciP2)(g2fc2)((fcip2)(pigi) + (pi(fci -P2))(gi(fci -P2))) 

(27r)8(((px +P2y - M^y + M^rl^Mq, + q^Y - Ml-Y + Ml,V^^) 

\qi\cPkid^k2d9 sint 



16(piP2)(((pi +P2 - - mg)2 + m2r2)ifcii|fc2i(|<fi| + m + |gi + Q2I cose) 

Here we also introduced the Tw, which are the widths of decay of W 
and V. Let us also note that |fci| and \k2 \ does not produce any divergence 
because if any of them is close to zero, then the appropriate scalar product 
(feiP2) or {q2k2) is also close to zero. 

The amplitude of the reaction represented in Figure [5] is equal to 



-/ X ig M (1 + 7^) / ^ -i ( (Pi +P2)m(Pi +P2) 

«(P2)— — ^ — "bOz— r-r« — T72- 9^- ^72 



V2' 2 '^''''{p,+p2Y-M^V"'" M, 



'^^^^^^^ 2 (pi+P2-fci)2-m2^^ 2 ^^^'^ 

(gi + g2)A(gi + g2)p ^ (l + 7!),(^^) (7) 



(gi + g2)2 - Ml- Ml, J ^ 

Its squared amphtude can be transformed to 

|,,r|2^ 8g^fclP2)(pigi)(g2fc2) 

' ' ((pi+P2)2-M^)2((Qi+Q2)2-M2,)((pi+p2-fci)^-m?) 

(This answer is also obtained by using mathematica) . In the same way 
one can show that the cross-section for this reaction is: 

, ^ 8/(fcip2)(pigi)(g2fc2) 

"'^ (27r)8(((pi +p2y - M^Y + M^Tl^mqr + 52)^ - M^^Y + Ml^V^^) 

\qi\d'^kid^k2d9 sin 9d(p 

16(pip2)(((pi +P2 - kiY - mlY + ■mlTl)\ki\\k2\{\qi\ + \q2\ + \qi + qi] cos 6*) 

3 The process as a decay chain. 

As we said in the introduction, instead of considering the whole process, 
we can assume that after the production in the collision of the quark and 
antiquark the boson decays into the lepton and heavy neutrino N and 
then the heavy neutrino N decays into the lepton and the quark-antiquark 
pair. The diagram in Figure [3] shows the last part of this process (the first 
part is trivial since we consider here the production of Wr at rest). 

In this section we calculate the width of this decay (it does not depend 
on the type of the neutrino). The expression for the amplitude is the 
following: 



V2 2 (gi+q2)2-M2, 

^„ (gi +g2)M(gi + 92),. ^ _ ig ^ (l + 7^) „. 



(6) 



(9) 
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Figure 3: The process represented as two consecutive decays 



Its squared modul can be simplified to 

|MP = 2g\Lq,){k.2q2) , . 

' ' {{{qi + q2Y-Ml,Y + M^Tl,y ^ > 

where L is the momentum of the neutrino and the other momenta are 
the same as in the previous section. We substitute this formula into the 
general expression for the width of decay ^ 

du, = ^J^\k, + q,+q,^L)\M\ ^g^ofcO^OgO (12) 

Again the delta-function of spatial momenta is taken off by integration 
over q2. The delta-function of energy is also taken off in the same way. 
We rewrite 

d'^ifi = \qx\^ d\q-i\ sin 6 d9 dip , 

where 6 and tp is polar and axial angle of q\ relative to the vector L — k2, 
and 



L°~k^^~q1-q^ = +M3-\k2\-\qi\-y \L - + \qi\^ - 2\q2\\L - fcaj cos( 

After integration over we obtain that 

_ \M\^\q-i_\sined-^k2dedp 
^ ^ (27r)n6(l(ji| + \q2\ + \qi + q2\ cos e)\k2 + qi +g2||A;2 
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4 Description of the programs 



We compare the angular distributions of the final state particles obtained 
using the full and simplified calculations. Two different programs per- 
form the Monte-Carlo simulation of the process using the von Neumann's 
method. The first program uses the cross-section formula obtained in sec- 
tion [5J the second one uses the simplified approach considered in section 
[3l More details of this simulation follow. 

In this estimation we assume that the masses of Wr and heavy neu- 
trino are equal to 1000 GeV and 500 GeV correspondingly. These are 
typical values for the searches of these new particles at the LHC. For 
the moment we assume that the initial quarks have the same energy and 
opposite momentum directions so that the Wr is produced at rest: 

pl^{E,0,0,E}, 
p2 = {E,0,0,-E}, (14) 

where E is the initial energy of quarks. 

The first program has two similar versions: the first one considers 
the case of Dirac-like neutrino where only the reaction 1 can take place, 
the considers considers the case of Majorana-like neutrino where both 
reactions 1 and 2 happen with a probability of 50% each. 

The maximum of the cross section is determined (maximization step) 
by the Monte Carlo method before the simulation step. At each iteration 
of the corresponding loops the final state particles 4-momenta take random 
values so that the energy and momentum are conserved. At the simulation 
step 2 or 4 loops are performed because antiquark can have positive or 
negative Pz and because reactions 1 and 2 are possible in the Majorana 
neutrino case. 

5 Results 

In Figure |4] one can see the results. They are represented as a comparison 
of the distributions of the minimal distance R in the rj — (j) plane from 
one of the final state leptons to one of the final state quarks. This dis- 
tribution is important because in real collider experiments the efficiency 
of reconstruction drops when this distance is small (to exact zero when 
it is smaller than 0.5). The shape is different. To estimate the system- 
atic error arising from the simulation as a decay chain we calculated the 
ratio of the number of events with J? > 1, where the drop of the lepton 
reconstruction efficiency due to the presence of a nearby jet is small, to 
the total number of events for the two simulations. They differ by 2.7%. 
This is our estimate of the systematic error of the simplified simulation. 
The difference between the cases of Dirac and Majorana heavy neutrinos 
is small (within 1%). 
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Eta-phi distance between a lepton and final quaric 
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Figure 4: Distributions of the minimal distance in the rj — (j) plane from one of 
the final state leptons to one of the final state quarks. Histogram - from the 
simulation as a decay chain (statistical errors arc negligible), points with errors 
- from the full matrix element with a Majorana heavy neutrino. 
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